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Recent progresses of the proton structure measurements and determination of parton distribution functions by ep 
collisions at HERA are introduced. 
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1. INTRODUCTION 

s : 

' Deep inelastic scattering (DIS) of leptons off nucleons has been the key for our understanding of the structure of 

j . the nucleon. HERA at DESY is a unique facility for colliding electrons (or positrons) with protons. Compared with 

£S| previous fixed-target DIS experiments, the large center-of-mass energy of about 320 GeV allows an extension of the 
exploitable kinematic phase space by two orders of magnitude both in Q 2 , the negative of the four-momentum transfer 
^ squared, and in the Bjorken scaling variable x. The maximum Q 2 reaches almost 10 5 GeV 2 , which corresponds to a 
. spatial resolution of 10~ 16 cm. The minimum x, which means the momentum fraction carried by the struck parton 1 , 
Qh) reaches almost 10 -5 . The cross section of the neutral-current (NC) DIS interaction, e ± p — > e ± A, can be written as 

i where a is the fine structure constant, Y± — 1 ± (1 — y) 2 with y being the inelasticity, and F%, xF$ and Fl are 
the structure functions of the proton. In the framework of the perturbative QCD inspired quark parton model, 
the structure functions can be directly related to the parton distribution functions (PDFs) which are probability 
densities of partons existing inside proton. At low Q 2 , dominantly contributing to the cross section is F2 which is an 
electric-charge squared weighted sum of all flavor quark PDFs. In the low- a; region, F% is dominated by sea-quark 
PDFs, and the DGLAP evolution of QCD ascribes the <3 2 -dependence of Fi ("scaling violation") as largely owing 
to gluon splitting into gg-pairs. Thus, HERA data provide crucial information on the small- a; sea-quark and gluon 
PDFs. At large Q 2 , xF% becomes significant, and gives information on valence quark PDFs. The structure function 
• i-H , Fl is zero in the naive quark-parton model, i.e. without QCD, but in leading order QCD, a finite value of Fl is 
expected in the small x region by being directly related to the gluon PDF. The cross section of the charged-current 
■ (CC) DIS interaction, e + (e~)p — > V(v)X, can be written as 

_ _ 1 

d 2 a(e + p) G% Mh, r/ . ,„ .,, , 
d 2 cr(e~p) Gp Mtr . ,„ , 9/ - 

= ^w^? [{u+c)+{l - v){d+1)l (3) 

thus bringing flavor sensitivity of the valence quark PDFs at large a;. 

Based on about 100 (20) pb" 1 of e+p (e~p) data collected until the year 2000 (HERA-I), both the HI and ZEUS 
experiments at HERA have measured the inclusive e ± p NC and CC DIS cross sections [1-10]. After the year 2000, 
HERA underwent a major upgrade (HERA-II) aiming for higher luminosity 2 , and until March 2007, HERA provided 
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1 In the proton's infinite momentum frame. 

2 The other aim of HERA-II is to provide collisions with longitudinally polarized beams, giving a direct sensitivity to the helicity 
structure of the Electro- Weak interaction. 



about 500 pb" 1 of e ± p collisions to each HI and ZEUS experiments. The increased statistics provides more sensitivity 
to the PDFs at large x and large Q 2 , in regions where the precisions of the HERA-I measurements are still statistically 
limited. In addition, the precision of the structure function measurements by tagging heavy-quarks in the final state 
can be significantly improved by making use of the increased statistics of HERA-II and also increased capability of 
heavy-quark identification with newly installed tracking devices. During March to June in the year 2007, HERA made 
a series of dedicated runs with reduced proton beam energies of 460 and 575 GeV as compared to the nominal one 
of 920 GeV. These data sets are essential in the first direct measurement of Fl, as will be discussed later. Although 
HERA ended its 15 years spanning operation history in June 2007, analyses using full data sets are lively ongoing, 
and significant progresses are being made. In this manuscript, recent results from HERA are presented, and their 
impacts are discussed. 

2. H1 AND ZEUS COMBINED ANALYSIS 
2.1. Combined cross sections 

PDFs are usually determined in global QCD analyses of DIS data, both from HERA and fixed-target experiments 
as well as jet production data from the TEVATRON [11, 12]. Given that high-Q 2 HERA data which are sensitive 
to large x PDFs are now available, HI and ZEUS used their own data alone to make PDF fits (called HI PDF 2000 
and ZEUS- JETS fits, respectively) [13, 14]. Statistical uncertainties of these data are very small at low Q 2 such 
that total uncertainties are dominated by systematic uncertainties, and data points are correlated through common 
systematic uncertainties, both within and across the data sets. Thus, consideration of point-to-point correlations 
between systematic uncertainties was essential. In fact, this issue is faced by all of such QCD analyses; there are 
basically two methods, the Hessian and Offset methods, but each method has its own advantages and shortcomings, 
and there is no agreed standard (see e.g. [15]). For example, in the Hessian method, each systematic error source is 
treated as an additional fit parameter which is optimized by the model provided by QCD. However, difficulties arise 
when there are inconsistencies between data sets, leading to shifts of systematic error parameters by far more than 
one standard deviation 3 . 

These drawbacks can be significantly reduced by averaging the data sets in a model-independent way prior to 
performing a QCD analysis on them. Recently, the HI and ZEUS experiments have made a combined analysis on 
their data in order to obtain such averaged data sets [16]. Almost all inclusive DIS cross section measurements with 
HERA-I data, including low Q 2 NC as well as high Q 2 NC and CC, are combined [1-10]. Since both sets of cross 
section data from HI and ZEUS measure the same true cross-sections, the combination procedure is a ^-minimization 
in which the parameters are the true values of the cross-section at each (x, Q 2 ) point and the correlated systematic 
error parameters. These parameters are thus determined in a model independent way using the Hessian method [17]. 
This combination procedure cross-calibrates the measurements, such that the overall uncertainty is reduced by much 
more than the simple average. In total, 1153 individual NC and CC measurements are averaged to 554 unique points, 
and 43 sources of correlated systematic uncertainty parameters are fitted. This yields a good quality of fit with the 
X 2 /cfo/=510/599. In the small Q 2 region of Q 2 < 12 GeV 2 the combined data set has precision much better than 
2%, whereas HI and ZEUS data each have a precision of ~ 3 %. Figure 1 illustrates the HI and ZEUS combined 
cross sections for three selected x bins as a function of Q 2 . The HI (open circles) and ZEUS data (open squared) 
are compared to the HI and ZEUS combined data (closed points). Measurements from the individual experiments 
have been displaced sideways for clarity. The error bars show the total uncertainty. At low Q 2 where the data are 
limited by systematic uncertainties, improvement in the total error is impressive. At higher Q 2 the combined data 
exhibit far smaller fluctuations as compared to individual data due to the increased statistical accuracy. 



3 Some global QCD analyses thus use non-statistical criteria to estimate the PDF uncertainties (Sx 2 = T 2 » 1). 
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Figure 1: HI and ZEUS combined cross sections for three selected x bins as a function Q 2 , together with the HI PDF 2000 
and ZEUS- JETS predictions. 



2.2. Combined PDF: HERAPDF 

A next-to-leading order (NLO) QCD analysis was performed using the combined HERA-I data set of e ± p NC and 
CC cross sections presented in the previous subsection as the sole input (called HERAPDF0.1) [18]. The PDFs are 
parameterized at Q\ = 4 GeV 2 with 11 free parameters. The constraints of the number sum-rules and the momentum 
sum rule are applied. All the data have values of the invariant mass of the hadronic system, W 2 > 300 GeV 2 and 
of x < 0.65 such that there is no need to account for target mass corrections or higher twist effects. Also, a cut of 
Q 2 > 3.5 GeV 2 is applied to avoid non perturbative effects. As the correlated systematic uncertainties are no longer 
crucial in the HERA-I combined data sets, different handling of systematic uncertainties resulted in similar ht results. 
For the central fit, the 43 systematic uncertainties from the individual HI and ZEUS data sets are quadratically 
combined to the statistical uncertainties, and the 4 sources of uncertainty from the combination procedure is treated 
with the Offset method. The x 2 P er degree of freedom of the fit was 477/562. 

Figure 2 shows the combined HERA-I cross sections together with the HERAPDF0.1 fit results. The fit gives a 
good description both to NC and CC cross sections, and both for collisions. Figure 3 shows the HERAPDF0.1 at 
Q 2 = 10 GeV 2 (right plot). Plotted together in the left plot are the HI 2000 PDF and ZEUS-JET PDFs, i.e. obtained 
by HI and ZEUS fitting to their own data. They have rather different shapes, and the size of the uncertainties 
are compatible to those of the global PDF analyses. This clearly shows that HERAPDF0.1 has a much better 
determination benefiting from the high precision of the combined data sets in which inconsistencies between the two 
data sets are resolved by the cross calibration. It is worth noting that the large amount of HERA-II data are not yet 
included in the combined data sets, and thus further improvements are foreseen for valence quark PDFs at large x, 
in particular for the d-quark PDF which will be determined from e + p CC data. Figure 4 presents the HERAPDF0.1 
at large Q 2 of 10000 GeV 2 compared to one of the global PDF analyses (CTEQ 6.5 [12] is chosen here), showing 
that the improvement in precision of PDFs holds until large Q 2 . The HERAPDF0.1 has a significant impact on the 
precision of the predictions for Standard Model processes such as W and Z production at LHC [19]. This due to the 
improved precision of the sea and the gluon PDFs at Q 2 ~ 10000 GeV 2 and 5 x 10 -4 <i<5x 10~ 2 . In fact, it has 
been shown that the precision of the predicted cross-sections at LHC is less than 2% in the mid-rapidity range for 
HERAPDF0.1 PDFs as compared to 15% for pre-HERA PDFs [20]. 
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Figure 3: HI PDF 2000 and ZEUS- JETS PDFs (left) and HERAPDF0.1 PDFs (right) at Q 2 = 10 GeV 2 . 



3. DIRECT MEASUREMENT OF F L 



Despite of the remarkable success of QCD as seen in previous sections, there are still several issues to be clarified. 
For example, validity of the DGLAP formalism should be confirmed more explicitly at small x. This is because the 
DGLAP equation is a leading log re-summation of In Q 2 terms, and thus it cannot be expected to work at small 
x where large ln(l/x) terms will emerge. Notice also that the gluon PDF is extracted 'indirectly' through the Q 2 - 
dependence of inclusive cross sections, and thus should be further extensively verified with other processes which 



HI and ZEUS Combined PDF Fit 




1 



10 J 



10" 



10' 



1 



X 



Figure 4: HERAPDF0.1 at Q 2 = 10000 GeV 2 compared to CTEQ6.5 PDF. 



have more direct sensitivity. In both of these senses, a measurement of Fl is of importance. As Fl is zero in a static 
view of proton, it is related to the dynamical picture of the proton. Therefore, Fl is in particular interesting at small 
x, and only HERA can perform such measurement. In the DGLAP formalism Fl is related to the gluon PDF at 
small x, thus giving a confirmation of indirectly determined gluon PDF in QCD analyses. 

The structure function Fl can be extracted from sets of cross sections measured with different center-of-mass 
energies, s, as follows: if we present cross sections in a reduced form ay defined as 



Fl can be extracted from measurements of ay made at various s, from the slope of a r with y 2 /Y+ 4 . Note that for 
measurements at fixed (x, Q 2 ) variation of y requires variation of s, since Q 2 = xys. The precision of the extracted 
Fl is dominated by the size of the 'lever arm' spanned in the axis of y 2 /Y+, and thus, extending the measurement 
up to y as high as possible is crucial. At low Q 2 , high y values correspond to low values of the scattered electron 
energy. As small energy depositions can also be caused by hadronic final state particles and hence fake electrons, 
there are large backgrounds expected at high y dominantly due to photoproduction precesses at Q 2 ~ GeV 2 . 
In HERA-II, HI and ZEUS has developed new measurement techniques optimized for large y. These were then 
validated with nominal beam energy data [21, 22] before starting the dedicated runs with reduced proton energies. 
In the HI analysis, photoproduction backgrounds are statistically subtracted by using the sample in which wrong 
sign of charge is reconstructed for the track linked to the electron candidates measured in calorimeter. A small 
correction to charge asymmetry in backgrounds was obtained from comparisons of samples of negative tracks in e + p 
scattering with samples of positive tracks in e~p scattering. In the ZEUS analysis, a part of the photoproduction 
backgrounds are tagged by using a special detector located close to the beampipe, and are used to understand and 
normalize the background MC samples. These experimental techniques are used in the analyses of special runs with 
reduced proton beam energies, i.e. in the Fl measurements. HI has measured Fl in middle to high Q 2 region, 
12 < Q 2 < 800 GeV 2 [23, 24]. Figure 5 shows the reduced cross section measured with the three different proton 
beam energies, plotted as function of y 2 /Y + . The data clearly shows a finite value of slope, hence Fl- Figure 6 
shows the extracted Fl as a function of x in various Q 2 bins. The result is consistent with the prediction obtained 
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4 Contribution from xF^ can be neglected at small Q 2 . 
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Figure 5: Reduced cross section oy at Q 2 — 25 GeV 2 as a function of y 2 /Y+, as measured by HI. 
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Figure 6: Ft as a function of x at various Q , as measured by HI. 



with the HI PDF 2000 fit. The values on Fl resulting from averages over x at fixed Q 2 are shown in Figure 7. Within 
the experimental uncertainties, the data are consistent with the QCD predictions, indicating the applicability of the 
DGLAP formalism at small x at HERA. While, Figure 8 shows the Fl as measured by ZEUS [25]. The measurement 
is consistent with the ZEUS-JETS prediction of Fl, i.e. a QCD prediction, but also with Fl = 0, within the current 
precision of the measurement. 



4. HEAVY QUARK CONTRIBUTION TO F 2 

In leading order of QCD, the heavy flavor production in ep collisions is dominated by the photon gluon fusion 
process (BGF), 75 — » QQ. Hence, it is sensitive to the gluon PDF and allows tests of the universality of the gluon 
PDF. The heavy quark contribution to F%, F% , can be obtained by extrapolating the visible measured cross section 
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Figure 7: Fl as a function of Q 2 at given values of x, compared with QCD predictions. 
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Figure 8: Fl as a function of x at various Q 2 , as measured by ZEUS. 



of heavy quark production to the full phase space using NLO QCD calculations. 

Figure 9 shows the charm contribution to F2, F^ , as measured by HI and ZEUS [26-34]. By using about 160 pb _1 
of luminosity of HERA-II data, ZEUS has measured the D*^ meson production cross section by tagging it with slow 
7r [26]. Production of Z? ± and D° mesons have also been measured making use of the capability of secondary vertex 
identification of the Micro Vertex Detector (MVD) which was installed at HERA-II [27]. A signed 2-dimensional 
decay length significance with respect to the beamspot measured in MVD gives a significant improvement in the 
signal to noise ratio, and thus in the statistical uncertainty of the measured cross section. These analyses are based 
on HERA-II data of about 130 pb _1 of luminosity. The results are plotted in the left plot of Figure 9, together with 
previous measurements by using HERA-I data, i.e. without the advantage of the MVD [28]. HI has measured D*^ 
meson production [29], as well as the inclusive charm and beauty production [30] cross sections, based on HERA-II 
data of about 350 and 190 pb _1 of luminosities, respectively. In the inclusive analysis, the charm and beauty fractions 
are determined using a neural network which includes the impact parameter of tracks to the primary vertex and the 
position of the secondary vertex as measured using the Vertex Detector (VTX) . These results are also shown in the 
left plot of Figure 9, compared to previous measurements using HERA-I data set [31-33]. The FSf measurement 
using the VTX gives a consistent result with that measured using the technique of tagging D*^ meson, thus providing 
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Figure 9: Heavy quark contribution to the structure function F2: charm contribution F| c (left) and bottom contribution F$ 
(right). 



a proof of the VTX analysis. 

In the right plot of Figure 9 shown is the bottom-quark contribution to F 2 , F^, as measured by HI in the inclusive 
analysis using VTX [30, 32, 33]. Also plotted is the ZEUS measurement tagging bottom quarks by requiring muon 
associated with jet, based on HERA-II data sample of about 40 pb" 1 [34]. In these F 2 C and F 2 b measurements, it 
is shown that a significant fraction of F 2 is owed by heavy quarks, about 30% by charm and about 3% by bottom, 
increasing with Q 2 . Precision of determining heavy-quark PDFs will be further improved by the analyses with full 
HERA-II data. 



5. MEASUREMENT OF xF 3 

The difference of cross sections between e + p and e~p collision data determines xF% as seen in Eq. 1. At HERA, 
xF% is, to a very good approximation, dominated only by the j-Z interference term, xF^ z , as 

Y 

cr r (e~) - o>(e + ) = -2— a e k z xF^ , (5) 
Y + 

where oy is the reduced cross section defined in Section 3, kz — . ■ -i a 1 — rs — 02T » 71 with Mz being the mass of Z 

4 sin "w cos "w leg -\- ivi g 

boson and 6\y being the Weinberg angle, and a e = —1/2 is the axial- vector coupling of electron to the Z boson. In 
leading order QCD, xF^ z can be written as 

xF^ z = 2x[e u a u uv + eddddv], (6) 

with an assumption that A u = (u sea — u + c — c) and = (d sea — d + s — s) are zero, where e u and a u are the charge 
and axial-vector couplings to the Z of the it-quark, and and ad are the corresponding quantities for the d-quark. 
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Figure 10: The f-Z interference term of the structure function xFz, xF^ z . 
A sum rule holds in leading order as 

f xF f<^ = 1 [\2uy+d v )dx=\. (7) 
JO x " Jo & 

The structure function xF^ z is thus determined by the valence quark PDFs and predicted to be only very weakly 
depending on Q 2 . 

Figure 10 shows the xF^ z measured by using both HI and ZEUS high Q 2 NC cross sections from HERA-I and 
HERA-II [35]. The measurement is well described by the HI 2000 PDF and ZEUS-JETS fits. The small difference 
between the HI and ZEUS fits indicates weak constraints on the behavior of the valence quark PDFs at smaller x. 
In the range of acceptance, the integral of F^ z is measured to be 

,•0.65 

/ F2 Z = 1.21 ± 0.09(stat) ± 0.08(syst), (8) 

JO. 02 

which is consistent with the results of the HI and ZEUS QCD fits, 1.12 ± 0.02 and 1.06 ± 0.02, respectively. 

6. SUMMARY 

HERA has provided the most precise inclusive structure function measurements, which has brought significant 
improvements to our knowledge on proton structure. For further comprehensive understanding, analyses with new 
techniques and with the full amount of data are being performed, and new high precision results are being produced. 
Final publications with ultimate precision will come in the next years. 
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